Neolithic genomic data from southern France showcase intensified interactions with hunter-gatherer communities.

Archaeological research shows that the dispersal of the Neolithic took a more complex turn when reaching western Europe, painting a contrasted picture of interactions between autochthonous hunter-gatherers (HGs) and incoming farmers. In order to clarify the mode, the intensity, and the regional variability of biological exchanges implied in these processes, we report new palaeogenomic data from Occitanie, a key region in Southern France. Genomic data from 28 individuals originating from six sites spanning from c. 5,500 to c. 2,500 BCE allow us to characterize regional patterns of ancestries throughout the Neolithic period. Results highlight major differences between the Mediterranean and Continental Neolithic expansion routes regarding both migration and interaction processes. High proportions of HG ancestry in both Early and Late Neolithic groups in Southern France support multiple pulses of inter-group gene flow throughout time and space and confirm the need for regional studies to address the complexity of the processes involved.

Origin and mobility of Iron Age Gaulish groups in present-day France revealed through archaeogenomics.

The Iron Age period occupies an important place in French history because the Gauls are regularly presented as the direct ancestors of the extant French population. We documented here the genomic diversity of Iron Age communities originating from six French regions. The 49 acquired genomes permitted us to highlight an absence of discontinuity between Bronze Age and Iron Age groups in France, lending support to a cultural transition linked to progressive local economic changes rather than to a massive influx of allochthone groups. Genomic analyses revealed strong genetic homogeneity among the regional groups associated with distinct archaeological cultures. This genomic homogenization appears to be linked to individuals' mobility between regions and gene flow with neighbouring groups from England and Spain. Thus, the results globally support a common genomic legacy for the Iron Age population of modern-day France that could be linked to recurrent gene flow between culturally differentiated communities.

Evidence for early dispersal of domestic sheep into Central Asia.

The development and dispersal of agropastoralism transformed the cultural and ecological landscapes of the Old World, but little is known about when or how this process first impacted Central Asia. Here, we present archaeological and biomolecular evidence from Obishir V in southern Kyrgyzstan, establishing the presence of domesticated sheep by ca. 6,000 BCE. Zooarchaeological and collagen peptide mass fingerprinting show exploitation of Ovis and Capra, while cementum analysis of intact teeth implicates possible pastoral slaughter during the fall season. Most significantly, ancient DNA reveals these directly dated specimens as the domestic O. aries, within the genetic diversity of domesticated sheep lineages. Together, these results provide the earliest evidence for the use of livestock in the mountains of the Ferghana Valley, predating previous evidence by 3,000 years and suggesting that domestic animal economies reached the mountains of interior Central Asia far earlier than previously recognized.

Ancient genomes from present-day France unveil 7,000 years of its demographic history.

Genomic studies conducted on ancient individuals across Europe have revealed how migrations have contributed to its present genetic landscape, but the territory of present-day France has yet to be connected to the broader European picture. We generated a large dataset comprising the complete mitochondrial genomes, Y-chromosome markers, and genotypes of a number of nuclear loci of interest of 243 individuals sampled across present-day France over a period spanning 7,000 y, complemented with a partially overlapping dataset of 58 low-coverage genomes. This panel provides a high-resolution transect of the dynamics of maternal and paternal lineages in France as well as of autosomal genotypes. Parental lineages and genomic data both revealed demographic patterns in France for the Neolithic and Bronze Age transitions consistent with neighboring regions, first with a migration wave of Anatolian farmers followed by varying degrees of admixture with autochthonous hunter-gatherers, and then substantial gene flow from individuals deriving part of their ancestry from the Pontic steppe at the onset of the Bronze Age. Our data have also highlighted the persistence of Magdalenian-associated ancestry in hunter-gatherer populations outside of Spain and thus provide arguments for an expansion of these populations at the end of the Paleolithic Period more northerly than what has been described so far. Finally, no major demographic changes were detected during the transition between the Bronze and Iron Ages.

Ancient genomic changes associated with domestication of the horse.

The genomic changes underlying both early and late stages of horse domestication remain largely unknown. We examined the genomes of 14 early domestic horses from the Bronze and Iron Ages, dating to between ~4.1 and 2.3 thousand years before present. We find early domestication selection patterns supporting the neural crest hypothesis, which provides a unified developmental origin for common domestic traits. Within the past 2.3 thousand years, horses lost genetic diversity and archaic DNA tracts introgressed from a now-extinct lineage. They accumulated deleterious mutations later than expected under the cost-of-domestication hypothesis, probably because of breeding from limited numbers of stallions. We also reveal that Iron Age Scythian steppe nomads implemented breeding strategies involving no detectable inbreeding and selection for coat-color variation and robust forelimbs.

Taming the late Quaternary phylogeography of the Eurasiatic wild ass through ancient and modern DNA.

Taxonomic over-splitting of extinct or endangered taxa, due to an incomplete knowledge of both skeletal morphological variability and the geographical ranges of past populations, continues to confuse the link between isolated extant populations and their ancestors. This is particularly problematic with the genus Equus. To more reliably determine the evolution and phylogeographic history of the endangered Asiatic wild ass, we studied the genetic diversity and inter-relationships of both extinct and extant populations over the last 100,000 years, including samples throughout its previous range from Western Europe to Southwest and East Asia. Using 229 bp of the mitochondrial hypervariable region, an approach which allowed the inclusion of information from extremely poorly preserved ancient samples, we classify all non-African wild asses into eleven clades that show a clear phylogeographic structure revealing their phylogenetic history. This study places the extinct European wild ass, E. hydruntinus, the phylogeny of which has been debated since the end of the 19th century, into its phylogenetic context within the Asiatic wild asses and reveals recent mitochondrial introgression between populations currently regarded as separate species. The phylogeographic organization of clades resulting from these efforts can be used not only to improve future taxonomic determination of a poorly characterized group of equids, but also to identify historic ranges, interbreeding events between various populations, and the impact of ancient climatic changes. In addition, appropriately placing extant relict populations into a broader phylogeographic and genetic context can better inform ongoing conservation strategies for this highly-endangered species.

Experimental conditions improving in-solution target enrichment for ancient DNA.

High-throughput sequencing has dramatically fostered ancient DNA research in recent years. Shotgun sequencing, however, does not necessarily appear as the best-suited approach due to the extensive contamination of samples with exogenous environmental microbial DNA. DNA capture-enrichment methods represent cost-effective alternatives that increase the sequencing focus on the endogenous fraction, whether it is from mitochondrial or nuclear genomes, or parts thereof. Here, we explored experimental parameters that could impact the efficacy of MYbaits in-solution capture assays of ~5000 nuclear loci or the whole genome. We found that varying quantities of the starting probes had only moderate effect on capture outcomes. Starting DNA, probe tiling, the hybridization temperature and the proportion of endogenous DNA all affected the assay, however. Additionally, probe features such as their GC content, number of CpG dinucleotides, sequence complexity and entropy and self-annealing properties need to be carefully addressed during the design stage of the capture assay. The experimental conditions and probe molecular features identified in this study will improve the recovery of genetic information extracted from degraded and ancient remains.

A New High-Throughput Approach to Genotype Ancient Human Gastrointestinal Parasites.

Human gastrointestinal parasites are good indicators for hygienic conditions and health status of past and present individuals and communities. While microscopic analysis of eggs in sediments of archeological sites often allows their taxonomic identification, this method is rarely effective at the species level, and requires both the survival of intact eggs and their proper identification. Genotyping via PCR-based approaches has the potential to achieve a precise species-level taxonomic determination. However, so far it has mostly been applied to individual eggs isolated from archeological samples. To increase the throughput and taxonomic accuracy, as well as reduce costs of genotyping methods, we adapted a PCR-based approach coupled with next-generation sequencing to perform precise taxonomic identification of parasitic helminths directly from archeological sediments. Our study of twenty-five 100 to 7,200 year-old archeological samples proved this to be a powerful, reliable and efficient approach for species determination even in the absence of preserved eggs, either as a stand-alone method or as a complement to microscopic studies.

Analysis of Ancient DNA in Microbial Ecology.

The development of next-generation sequencing has led to a breakthrough in the analysis of ancient genomes, and the subsequent genomic analyses of the skeletal remains of ancient humans have revolutionized the knowledge of the evolution of our species, including the discovery of a new hominin, and demonstrated admixtures with more distantly related archaic populations such as Neandertals and Denisovans. Moreover, it has also yielded novel insights into the evolution of ancient pathogens. The analysis of ancient microbial genomes allows the study of their recent evolution, presently over the last several millennia. These spectacular results have been attained despite the degradation of DNA after the death of the host, which results in very short DNA molecules that become increasingly damaged, only low quantities of which remain. The low quantity of ancient DNA molecules renders their analysis difficult and prone to contamination with modern DNA molecules, in particular via contamination from the reagents used in DNA purification and downstream analysis steps. Finally, the rare ancient molecules are diluted in environmental DNA originating from the soil microorganisms that colonize bones and teeth. Thus, ancient skeletal remains can share DNA profiles with environmental samples and identifying ancient microbial genomes among the more recent, presently poorly characterized, environmental microbiome is particularly challenging. Here, we describe the methods developed and/or in use in our laboratory to produce reliable and reproducible paleogenomic results from ancient skeletal remains that can be used to identify the presence of ancient microbiota.

Twenty-five thousand years of fluctuating selection on leopard complex spotting and congenital night blindness in horses.

Leopard complex spotting is inherited by the incompletely dominant locus, LP, which also causes congenital stationary night blindness in homozygous horses. We investigated an associated single nucleotide polymorphism in the TRPM1 gene in 96 archaeological bones from 31 localities from Late Pleistocene (approx. 17 000 YBP) to medieval times. The first genetic evidence of LP spotting in Europe dates back to the Pleistocene. We tested for temporal changes in the LP associated allele frequency and estimated coefficients of selection by means of approximate Bayesian computation analyses. Our results show that at least some of the observed frequency changes are congruent with shifts in artificial selection pressure for the leopard complex spotting phenotype. In early domestic horses from Kirklareli-Kanligecit (Turkey) dating to 2700-2200 BC, a remarkably high number of leopard spotted horses (six of 10 individuals) was detected including one adult homozygote. However, LP seems to have largely disappeared during the late Bronze Age, suggesting selection against this phenotype in early domestic horses. During the Iron Age, LP reappeared, probably by reintroduction into the domestic gene pool from wild animals. This picture of alternating selective regimes might explain how genetic diversity was maintained in domestic animals despite selection for specific traits at different times.

Evidence for a retroviral insertion in TRPM1 as the cause of congenital stationary night blindness and leopard complex spotting in the horse.

Leopard complex spotting is a group of white spotting patterns in horses caused by an incompletely dominant gene (LP) where homozygotes (LP/LP) are also affected with congenital stationary night blindness. Previous studies implicated Transient Receptor Potential Cation Channel, Subfamily M, Member 1 (TRPM1) as the best candidate gene for both CSNB and LP. RNA-Seq data pinpointed a 1378 bp insertion in intron 1 of TRPM1 as the potential cause. This insertion, a long terminal repeat (LTR) of an endogenous retrovirus, was completely associated with LP, testing 511 horses (χ(2)=1022.00, p<<0.0005), and CSNB, testing 43 horses (χ(2)=43, p<<0.0005). The LTR was shown to disrupt TRPM1 transcription by premature poly-adenylation. Furthermore, while deleterious transposable element insertions should be quickly selected against the identification of this insertion in three ancient DNA samples suggests it has been maintained in the horse gene pool for at least 17,000 years. This study represents the first description of an LTR insertion being associated with both a pigmentation phenotype and an eye disorder.

From genes to phenotypes - evaluation of two methods for the SNP analysis in archaeological remains: pyrosequencing and competitive allele specific PCR (KASPar).

The amplification length of the DNA fragments is one major limitation of most paleogenetic analyses. Routinely, only fragments below 200 bp can be amplified, significantly reducing the content of genetic information. Although overlapping PCR strategies and next generation sequencing techniques have strongly improved data mining recently, these methods are still expensive and time consuming. In contrast, SNP analyses are easy to handle, fast and cheap. In this study, we compare two methods of SNP detection as to efficiency, cost and reliability for their use in ancient DNA applications: pyrosequencing and competitive allele specific PCR (KASPar). Our sample set consisted of 16 horse bones from two Scythian graves (600-800 BC). In conclusion, both approaches produced reliable results for most allelic patterns. But an indel of 11 bp (ASIP) could not be detected in the KASPar approach and produced problems in the pyrosequencing method (70% success rate). In such cases, we recommend checking allelic distribution using a gel approach or capillary sequencing. Overall, in comparison with the traditional mode of ancient DNA investigations (PCR, cloning, capillary sequencing), both approaches are superior for SNP analyses especially of large sample sets.

Genotypes of predomestic horses match phenotypes painted in Paleolithic works of cave art.

Archaeologists often argue whether Paleolithic works of art, cave paintings in particular, constitute reflections of the natural environment of humans at the time. They also debate the extent to which these paintings actually contain creative artistic expression, reflect the phenotypic variation of the surrounding environment, or focus on rare phenotypes. The famous paintings "The Dappled Horses of Pech-Merle," depicting spotted horses on the walls of a cave in Pech-Merle, France, date back ~25,000 y, but the coat pattern portrayed in these paintings is remarkably similar to a pattern known as "leopard" in modern horses. We have genotyped nine coat-color loci in 31 predomestic horses from Siberia, Eastern and Western Europe, and the Iberian Peninsula. Eighteen horses had bay coat color, seven were black, and six shared an allele associated with the leopard complex spotting (LP), representing the only spotted phenotype that has been discovered in wild, predomestic horses thus far. LP was detected in four Pleistocene and two Copper Age samples from Western and Eastern Europe, respectively. In contrast, this phenotype was absent from predomestic Siberian horses. Thus, all horse color phenotypes that seem to be distinguishable in cave paintings have now been found to exist in prehistoric horse populations, suggesting that cave paintings of this species represent remarkably realistic depictions of the animals shown. This finding lends support to hypotheses arguing that cave paintings might have contained less of a symbolic or transcendental connotation than often assumed.

Origin and history of mitochondrial DNA lineages in domestic horses.

Domestic horses represent a genetic paradox: although they have the greatest number of maternal lineages (mtDNA) of all domestic species, their paternal lineages are extremely homogeneous on the Y-chromosome. In order to address their huge mtDNA variation and the origin and history of maternal lineages in domestic horses, we analyzed 1961 partial d-loop sequences from 207 ancient remains and 1754 modern horses. The sample set ranged from Alaska and North East Siberia to the Iberian Peninsula and from the Late Pleistocene to modern times. We found a panmictic Late Pleistocene horse population ranging from Alaska to the Pyrenees. Later, during the Early Holocene and the Copper Age, more or less separated sub-populations are indicated for the Eurasian steppe region and Iberia. Our data suggest multiple domestications and introgressions of females especially during the Iron Age. Although all Eurasian regions contributed to the genetic pedigree of modern breeds, most haplotypes had their roots in Eastern Europe and Siberia. We found 87 ancient haplotypes (Pleistocene to Mediaeval Times); 56 of these haplotypes were also observed in domestic horses, although thus far only 39 haplotypes have been confirmed to survive in modern breeds. Thus, at least seventeen haplotypes of early domestic horses have become extinct during the last 5,500 years. It is concluded that the large diversity of mtDNA lineages is not a product of animal breeding but, in fact, represents ancestral variability.

An efficient multistrategy DNA decontamination procedure of PCR reagents for hypersensitive PCR applications.

BACKGROUND: PCR amplification of minute quantities of degraded DNA for ancient DNA research, forensic analyses, wildlife studies and ultrasensitive diagnostics is often hampered by contamination problems. The extent of these problems is inversely related to DNA concentration and target fragment size and concern (i) sample contamination, (ii) laboratory surface contamination, (iii) carry-over contamination, and (iv) contamination of reagents. METHODOLOGY/PRINCIPAL FINDINGS: Here we performed a quantitative evaluation of current decontamination methods for these last three sources of contamination, and developed a new procedure to eliminate contaminating DNA contained in PCR reagents. We observed that most current decontamination methods are either not efficient enough to degrade short contaminating DNA molecules, rendered inefficient by the reagents themselves, or interfere with the PCR when used at doses high enough to eliminate these molecules. We also show that efficient reagent decontamination can be achieved by using a combination of treatments adapted to different reagent categories. Our procedure involves γ- and UV-irradiation and treatment with a mutant recombinant heat-labile double-strand specific DNase from the Antarctic shrimp Pandalus borealis. Optimal performance of these treatments is achieved in narrow experimental conditions that have been precisely analyzed and defined herein. CONCLUSIONS/SIGNIFICANCE: There is not a single decontamination method valid for all possible contamination sources occurring in PCR reagents and in the molecular biology laboratory and most common decontamination methods are not efficient enough to decontaminate short DNA fragments of low concentration. We developed a versatile multistrategy decontamination procedure for PCR reagents. We demonstrate that this procedure allows efficient reagent decontamination while preserving the efficiency of PCR amplification of minute quantities of DNA.

Coat color variation at the beginning of horse domestication.

The transformation of wild animals into domestic ones available for human nutrition was a key prerequisite for modern human societies. However, no other domestic species has had such a substantial impact on the warfare, transportation, and communication capabilities of human societies as the horse. Here, we show that the analysis of ancient DNA targeting nuclear genes responsible for coat coloration allows us to shed light on the timing and place of horse domestication. We conclude that it is unlikely that horse domestication substantially predates the occurrence of coat color variation, which was found to begin around the third millennium before the common era.

Freshly excavated fossil bones are best for amplification of ancient DNA.

Despite the enormous potential of analyses of ancient DNA for phylogeographic studies of past populations, the impact these analyses, most of which are performed with fossil samples from natural history museum collections, has been limited to some extent by the inefficient recovery of ancient genetic material. Here we show that the standard storage conditions and/or treatments of fossil bones in these collections can be detrimental to DNA survival. Using a quantitative paleogenetic analysis of 247 herbivore fossil bones up to 50,000 years old and originating from 60 different archeological and paleontological contexts, we demonstrate that freshly excavated and nontreated unwashed bones contain six times more DNA and yield twice as many authentic DNA sequences as bones treated with standard procedures. This effect was even more pronounced with bones from one Neolithic site, where only freshly excavated bones yielded results. Finally, we compared the DNA content in the fossil bones of one animal, a approximately 3,200-year-old aurochs, excavated in two separate seasons 57 years apart. Whereas the washed museum-stored fossil bones did not permit any DNA amplification, all recently excavated bones yielded authentic aurochs sequences. We established that during the 57 years when the aurochs bones were stored in a collection, at least as much amplifiable DNA was lost as during the previous 3,200 years of burial. This result calls for a revision of the postexcavation treatment of fossil bones to better preserve the genetic heritage of past life forms.

Minimizing DNA contamination by using UNG-coupled quantitative real-time PCR on degraded DNA samples: application to ancient DNA studies.

PCR analyses of ancient and degraded DNA suffer from their extreme sensitivity to contamination by modern DNA originating, in particular, from carryover contamination with previously amplified or cloned material. Any strategy for limiting carryover contamination would also have to be compatible with the particular requirements of ancient DNA analyses. These include the need (i) to amplify short PCR products due to template fragmentation; (ii) to clone PCR products in order to track possible base misincorporation resulting from damaged templates; and (iii) to avoid incomplete decontamination causing artifactual sequence transformation. Here we show that the enzymatic decontamination procedures based upon dUTP- and uracil-N-glycosylase (UNG) can be adapted to meet the specific requirements of ancient DNA research. Thus, efficiency can be improved to vastly reduce the amplification of fragments < or = 100 bp. Secondly, the use of an Escherichia coli strain deficient in both UNG and dUTPase allows for the cloning of uracil-containing PCR products and offers protection from plasmid DNA contamination, and, lastly, PCR products amplified from UNG-degraded material are free of misleading sequence modifications.